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ABSTRACT 

We report results from a deep high-frequency search for pulsars within the central parsec of Sgr A* using the 
Green Bank Telescope. The observing frequency of 15 GHz was chosen to maximize the likelihood of detecting 
normal pulsars (i.e. with periods of ~ 500 ms and spectral indices of ~ — 1.7) close to Sgr A*, that might be used 
as probes of gravity in the strong-field regime; this is the highest frequency used for such pulsar searches of the 
Galactic Center to date. No convincing candidate was detected in the survey, with a 10er detection threshold of 
~ 10/.dy achieved in two separate observing sessions. This survey represents a significant improvement over 
previous searches for pulsars at the Galactic Center and would have detected a significant fraction (> 5%) of the 
pulsars around Sgr A*, if they had properties similar to those of the known population. Using our best current 
knowledge of the properties of the Galactic pulsar population and the scattering material toward Sgr A*, we 
estimate an upper limit of 90 normal pulsars in orbit within the central parsec of Sgr A*. 

Subject headings: Galaxy: center- pulsars: general 



1, INTRODUCTION 

The detection of radio-emitting neutron stars within the 
central parsec of our Galaxy would provide us with an un- 
precedented opportunity to study the super-massive black 
hole Sgr A* and its environs. For example, a single orbit- 
ing pulsar would yield a direct probe of the magneto-ionized 
accretion environment around a black hole, through mea- 
surements of t emporal changes in t he dispersion and rota- 
tion measures (Cord es & Lazioll997l) . Pulsars orbiting within 
the curved space-time around Sgr A* (with orbital periods of 
< 100 years) could serve as probes of gravity in the strong- 
field regime, at field strengths far larger than those accessible 
with neutron star binaries. The long-term timing of such pul- 
sars, supplemented by accurate astrometry, would allow pre- 
cise determination of their three-dimensional orbital motion 
around Sgr A*. Depending on the properties of the pulsars 
and their orbits, it should be possible to measure subtle gen- 
eral relativistic deviations from Keplerian orbits (e.g. time 
dilation, gravita tional redshifts, frame dragging, Shap iro de- 
lays, etc; e.g. ICordes et al.1l2004t iPfahl & Loebll20o4 . and it 
may even be possible to determine the spin of th e black hole 
(e.g. iWex & Kopeikhil999T: Kramer et al.ll2004l) . 

While theoretical estimates indicate that 100- 1000 radio 
pulsars with periods < 100 years should be orbiting Sgr A* 
(IPfahl & Loebl [2004), the observational evidence for neu- 
tron stars at the Galactic Center (GC) is mostly indirect. 
For example, recent studies have found a number of dense 
clusters of young, massive stars within ~ 1 pc of Sgr A * 
dSchodel et al.ll2003t iGhez et alj|2005t iPaumard et ai1l2006l) . 
while IWang et all ((2006) report X-ray observations of a pul- 
sar wind nebula near the massive stellar complex IRS 13, with 
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properties consistent with it being powered by a young neu- 
tron star. Long term monitoring by Chandra has revealed 
an excess of t ransient sou r ces w ithin a parsec of Sgr A*, 
interpreted by Mun o et all (120051) as a population of X-ray 
binaries. The flar ing radio and X-ray source detected by 
iBower et al.l (120051) . ~ 0.1 pc from Sgr A*, is also likely to 
be an X-ray binary. 

Despite the above evidence for massive stars around 
Sgr A*, there is a remarkable dearth of r adio pulsar detections 
there, despite several deep searches (e.g. IJohnston et al.l[l995T: 
Johnst on et al.ll2006t iDeneva et al.ll2009l) . The closest known 
radio pulsars are 11' from Sgr A*, and less than one percent of 
the known pulsar population lies within a degree of the Galac- 
tic Cen ter, despite indicati ons of a large population in its en- 
virons (IDeneva et al.ll2009l) . The reason for this deficit is well 
understood: hyper-strong scattering of radio waves by the tur- 
bulent, ionized gas within the central 100 pc of Sgr A*, which 
results in temporal smearing of pulsed signals. This pulse 
broadening has a strong frequency dependence, oc v~ A , mak- 
ing it near-impossible to detect pulsars at the typical ob serving 
frequencies of < 1.4 GHz (e.g. lLazio & Cord es 1998). 

To overcome the effects of temporal smearing, searches 
for pulsars at the GC have been carried out at progressively 
higher observing frequencies over the last few years (e.g. 
IJohnston etai] 12006; De neva et aI1l2009l) . albeit as yet with- 
out a detection in the central 25 pc. In this work, we report 
results from a deep Green Bank Telescope (GBT) search for 
pulsars toward Sgr A* at ~ 15 GHz, the highest observing 
frequency used till date. The choice of this frequency is moti- 
vated in Section $2] and the observations and results described 
in Section §|3] Finally, Section ^4] discusses the constraints 
placed by our observations on the GC pulsar population, and 
the prospects for pulsar detections in future surveys. 

2. A PULSAR SEARCH STRATEGY FOR THE GALACTIC CENTER 

Previous surveys of the GC pulsar population have concen- 
trated on maximizing the likel ihood of detection within a few 
degre es of Sgr A* (e.g. IJohnston et al.l 120061: IDeneva et al.l 
2009). This has the effect of lowering the optimal observing 
frequency, because a larger volume can be searched in a given 
amount of integration time at lower frequencies, due to the 
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larger telescope beam. However, the region containing pul- 
sars capable of probing gravitational effects in the strong-field 
regime is quite small, < 0.4", much smaller than the GBT 
beam at even frequencies ^ 10 GHz. Since our long-term 
goal is to use pulsars as probes of strong-field gravity around 
Sgr A*, we optimized our search strategy to detect pulsars in 
the immediate v i cinit y of (i.e. in orbit around) Sgr A* itself. 
ICordes & Laziol d 19971) have previously described the obser- 
vational challenges in detecting periodic emission at the GC; 
their main arguments are summarized below. 

The optimal observing frequency depends on three com- 
peting effects, the temporal scattering, the pulsar emission 
spectrum, and the telescope sensitivity, all with very differ- 
ent frequency dependencies. For example, the generally steep 
spectral energy dependence of the pulsar emission (S„ oc v a , 
with the distribution centered on a ~ -1 .7) drives searches to 
low frequencies, but competes against the strong v~ A spectral 
dependence of temporal smearing, due to which pulsars are 
severely attenuated at frequencies below which the scattering 
time exceeds the pulse period. We note, in passing, that the 
intervening ionized plasma also introduces a dispersive delay 
in the pulsar arrival times which can lead to a smearing of 
the pulse in individual filterbank frequency channels used to 
detect the signal, thus degrading pulse detectability. While 
this effect has been formally included in the present analy- 
sis, it is insignificant for the fine spectral resolution offered 
by modern filterbanks, for the estimated dispersion measure 
of ~ 1500- 3000 pc cm -3 for the GC region dLazio & Cordesl 
1998) 

The primary uncertainty in the temporal smearing time lies 
in the distance of the scattering screen from t he Galactic Cen- 
ter. Th e best estimate of this distance is from iLazio & Cordesl 
( 1998), who combined all known tracers of ionized gas (the 
scattering diameters of masers and OH/IR stars, free-free 
emission and absorption, etc.) in a maximum likelihood anal- 
ysis to obtain Z) scat = 133^0°pc. Making the usual simplify- 
ing assumption that that the scattering material is confined to 
a thin screen results in a temporal smearing timescale for pul- 
sars near Sgr A* of 

r ISM = 0.116 (ibzH-Y i^-Y sec. (1) 
VIOOpc / ^lOGHz/ V ' 

At a given frequency, a pulsar with period smaller than 
~ T scat would effectively have its pulses smeared into one an- 
other. In other words, for a given observing frequency, it is not 
possible to detect a pulsar whose period is much shorter than 
the temporal broadening time at that frequency (see Fig. [TJ. 
For instance, a pulsar of period 50 ms is dominated by tem- 
poral smearing at frequencies < 1 1 GHz, while the detection 
of a pulsar of period 5 ms is strongly hampered at frequencies 
< 26 GHz. Combining these effects, we estimate the tempo- 
ral smearing time for a pulse of intrinsic width Tj ntr j n sic, and 
with contributions from scattering (tismX interstellar disper- 
sion (tdm) and temporal binning (r res ) to the overall observed 
pulse width to be 

lobs = [intrinsic + T ISM + T DM + T res] • (2) 

This equation affords a sufficiently good approximation to 
the effects of temporal broadening, shown in detail in Fig.Q] 
when neither the intrinsic pulse profile nor the temporal 
broadening kernel are known in detail. 
Next, the steepness of the pulsar emission spectrum ad- 



versely impacts searches at high frequencies. The trend from 
younger to older pu lsars is a steepenin g of the spectral in- 
dex, from -1 to -2 (lLorimer et al.ll 19951) with a mea n of -1.6. 
Using a sample of 266 pulsars, iMaron et al.1 (120001) find that 
the flux density typically has a single power la w {S„ oc v a ), 
with a n average spectral index (a) = —1.8 ±0.2 (IMaron et al.l 
2000). For this paper we use a mean pulsar spectral in- 
dex of —1.7, a valu e int ermediate be t ween that derived by 
lLorimer et ail d 19951) and IMaron etail (12000). About 10% of 
the sample require a dual power law, with a steeper frequency 
dependence above ~ 1 .5 GHz; these pulsars would be difficult 
to detect at high frequencies. Conversely, there exists a minor- 
ity population (<2%) of pulsars that show a flat tening or even 
an upturn in their spectra at higher frequencies dKramer et al.l 
[19961 !LohmerinLll2008l) . 

Finally, the frequency dependencies of the performance of 
the telescope and its receivers, the Galactic background, and 
the sky brightness, all have an impact on the choice of ob- 
serving frequency. In the case of the GBT, the receiver con- 
tributions to the system temperature are roughly flat from 
1.5-20 GHz. At lower frequencies (v < 10 GHz), strong 
emission from the GC region is a significant (and sometimes 
overlooked) contributor t o the syste m temperature, with Tbg= 
340(i//2.7GHz)" 2 7 K dReich et al.ll 19901) . Conversely, the at- 
mosphere makes a significant contribution to the system tem- 
perature at high frequencies, > 15 GHz, especially at the low 
GC elevations when observed with the GBT. The lowest GBT 
system temperatures can be shown to arise in the frequency 
range ~ 10-18 GHz, with the upper and lower bounds deter- 
mined by contributions to T sys from the atmosphere and the 
GC region, respectively. 

All of the foregoing effects were combined to determine 
the optimal observing frequency for the GBT search for pul- 
sars at the Galactic Center. This was done by computing 
the signal-to-noise ratio (S/N) for a 10-hour GBT integration 
as a f unction of frequency (using equation 1 of iDewey et al.l 
119851) for a pulsar with a mean flux density of = lmJy 
at a frequency of 1 GHz, and an intrinsic pulsar duty cy- 
cle of 10%. This was done for three representative spectral 
indices, a = (—1.0 ,—1.7 ,-2.5) (where S v oc v a ). Finally, 
three representative periods (5, 50, and 500 ms) were also 
used, to illustrate the GBT's sensitivity to different pulsar 
populations (corresponding, roughly, to millisecond pulsars, 
partially-recycled or young pulsars, and normal pulsars, re- 
spectively). 

The results of this analysis are summarized in Fig. [2] where 
the S/N for a 10-hour GBT integration is plotted versus fre- 
quency; we use S/N = 10 as the detection threshold. It is clear 
from the figure that, for the assumed flux density (1 mJy at 
1 GHz), recycled or millisecond pulsars would not be de- 
tectable with the GBT unless they have flat spectral indices 
(a ~ -1), while pulsars with a = -2.5 would be entirely un- 
detectable. It is also clear that the optimal search frequency 
for "normal" pulsars (period ~ 500 ms, and a ~ —1.7) lies 
in the range 10-16 GHz, with the upper end of this fre- 
quency range also allowing the detection of shorter-period, 
flat-spectrum pulsars. The frequency range 12- 12.8 GHz is 
also affected by strong satellite-based interference at the GBT. 
This motivated our choice of 15 GHz as the optimal observ- 
ing frequency for the GBT search for pulsars at the Galactic 
Center. 

3. OBSERVATIONS, DATA ANALYSIS, AND RESULTS 
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FIG. 1 . — The effect of temporal broadening on a pulse with an intrinsic pulse width of 5% of the spin period, P S p m , assuming a one-sided exponential filter for 
the temporal broadening. The blue line indicates the fraction of the total pulsed flux density from the pulsar that remains pulsed after the scattering. The red line 
depicts the "Fourier SNR" (i.e. the sum of average noise-level normalized and mean-subtracted Fourier amplitudes) of a 16-harmonic summation of the power in 
a standard Fourier search. The latter indicates the detectability of pulsars in a blind survey and implicitly incorporates the pulse profile width increases and pulse 
shape changes caused by scattering. Most sensitivity is lost beyond T scat ~ 0.2P sp j„. 



TABLE 1 
GBT Observational Details 





fobs 


On-Src 




Sensitivity 


Date 


(GHz) 


Time (hrs) 


(K) a 


(/^y) b 


22 Jun 2006 


14.8 


2.5 


32 


18 


28 Jun 2006 


14.8 


4.5 


38 


16 


29 Jun 2006 


14.8 


5.5 


32 


12 


10 Aug 2006 


14.8 


3.0 


37 


19 


1 1 Aug 2006 


14.8 


2.25 


50 


30 


30 Aug 2008 


14.4 


4.25 


36 


16 


31 Aug 2008 


14.4 


5.5 


33 


13 



a The quoted T sys values are the averages for each run, af- 
ter correcting for the different elevations of the GC and the 
calibrators. 

b The last column lists the 10cr detection threshold flux den- 
sity for pulsars with a duty cycle of 10%, in each individual 
observing session. Note that the final searches combined 
data from adjacent days, and hence had a lower detection 
threshold flux density, ~ 10/^Jy apiece for the runs of 28 
and 29 June, 2006, and 30 and 31 August, 2008. 



3.1. The observations 

The search for pulsars at the Galactic Center was carried 
out with the Ku-band receiver of the GBT in the summers 
of 20 06 and 2008, with the GBT Pulsar Spigot (Kaplan et al.l 
12005b used as the backend. The observations were typically 



restricted to GC elevations higher than 10 degrees, to reduce 
atmospheric contributions to the system temperature (T sys ). 
The GC is only visible with the GBT for ~ 6 hours above this 
elevation limit per day, and the observations were therefore 
broken up into multiple short sessions. While attempts were 
made to group these sessions as close together as possible (so 
as to combine multiple sessions in a search), this was often 
not possible due to poor high-frequency observing weather. 

The 2006 observing runs were on June 22, 28 and 29, 
July 1 0, and August 1 and 1 1 . The July 1 data were affected 
by strong broad-band radio frequency interference (RFI) due 
to an overloaded network switch, and will not be discussed 
further; data from the other runs were not affected by this 
issue. The 2008 observations were on August 30 and 31, 
with each run preceded by a 5-hour observation of a blank 
field (J2000 co-ordinates: Right Ascension=llh46m08.1s, 
Declination=-27d47'32.9"), with the same observing param- 
eters as for the GC; the latter served to test for systematic 
effects in the data. The on-source times for the GC are listed 
in column (3) of TableQ] 

All observing runs included short observations of strong 
known pulsars (B 1800-21 in most of the 2006 runs and the 
magnetar, XTEJ18 10-197, in August 2006 and 2008), to test 
that the system was working properly; these were clearly de- 
tected in every run. 

All observing sessions used the maximum allowed Spigot 
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FIG. 2. — The frequency dependence of the detectability of pulsars at the Galactic Center for three representative periods, P=5ms (green), 50ms (red) and 
500 ms (blue) for a range of power-law spectral indices a=-\ (thin line), —1.7 (solid line), -2.5 (dashed line). The vertical axis gives the expected signal-to-noise 
ratio for a 10 hour integration with the GBT at the Galactic Center, normalized to a pulsar of period 1 s, flux density S u = 1 mjy at a frequency of 1 GHz, and 
whose spectrum has the frequency dependence S„ oc u a . The frequency range for the current suite of GBT receivers is shown at the top of the figure. 



bandwidth of 800 MHz, sub-divided into 1024 channels, with 
two circular polarizations, 3-level sampling, and a dump rate 
of 81.92 /its. The 800 MHz band was centered at 14.8 GHz in 
the 2006 runs, and at 14.4 GHz in the 2008 runs; the shift in 
the central frequency was to exclude some RFI (at 14.88 GHz) 
from the observing band. In 2006, an automatic level con- 
troller was used to hold the power levels fixed on the GC. 
Unfortunately, this was found to itself result in intermittent 
power jumps, and we hence chose to not use this feature in 
2008. In all runs, observations of strong nearby calibrators 
(or the GC itself) were used to correct the telescope pointing 
and surface, typically every 2-3 hours. System temperatures 
were measured on the calibrators by firing a noise diode, and 
ranged between 26 K and 40 K on most observing runs; how- 
ever, the short session on August 1 1 had T sys ~ 50 K on the 
calibrator, due to overcast conditions. After taking into ac- 



count the slightly different elevations of the GC and the cal- 
ibrators, the estimated average system temperatures towards 
the GC were 32-38 K (again, except on 11 August, where 
T sys ~ 50 K). The gain of the Ku-band receiver is 1.5 KJy" 1 . 
Finally, the size of the GBT beam at 14.8 GHz is ~ 50", cor- 
responding to a spatial radius of ~ 1 pc at the distance of the 
GC (7.9!°;* kpc: IReid et all2009l) . 

TableQjpresents a summary of the observational details; the 
columns are (1) the observing date, (2) the central frequency, 
(3) the GC on-source time, (4) the typical system temperature 
on the GC during the run, and (5) the 10er detection threshold 
for pulsars w ith a duty cycle of 10% in each session (using 
equation 1 of lDewey et al.l[l985T) . As discussed below, the fi- 
nal searches were carried out on the combined datasets from 
(1) June 28 and 29, 2006 (~ 10 hours), (2) August 10 and 
11, 2006 (~ 5.25 hours), and (3) August 30 and 31, 2008 
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(~ 9.75 hours). The first and third of these had the best sen- 
sitivity, with lOcr detection thresholds of ~ 10/iJy each, again 
assuming a pulsar duty cycle of 10%. 

3.2. Data analysis 

The observations were re duced and ana lyzed with the 
PRESTC0 software package dRansoml 1200 ll) . Initially, all 
data were time-domain clipped at the 6a level to remove any 
strong RFI. As noted above, the 2006 data suffered from ran- 
dom jumps in the power level due to an automatic level con- 
troller in the signal path. The effect of these jumps was re- 
moved by forcing data between jumps to a common mean of 
zero, with data in the immediate vicinity of the jump points 
edited out. Data from successive days (28 and 29 June 2006, 
10 and 1 1 August 2006, and 30 and 31 Aug 2008) were com- 
bined into a single data stream in order to enhance the de- 
tectability of faint pulsars during the periodicity search. The 
data from 28 and 29 June 2006 and 10 and 11 August 2006 
were also analysed independently, to examine the possibility 
of intermittent signals. The data streams were de-dispersed 
for dispersion measures of 1500, 2000, 2500, 3000, 3500, 
4000 and 5000 cm" 3 pc. Finer sampling in dispersion mea- 
sure is not required for the detection of slow pulsars be- 
cause the relative dispersion delay across the 800 MHz band 
at 14.8 GHz is only 2 ms for every 1000 cm" 3 pc. Note that 
our search was only sensitive to pulsars with periods < 2 s, 
due to power fluctuations associated with atmospheric vari- 
ability that reduced our sensitivity and precluded a search for 
pulsed emission at longer periods. 

The search for pulsars was carried out using standard tech- 
niques: the signal was Fourier transformed, a red-noise re- 
duction was applie d to the spectrum, an d the PRESTO algo- 
rithm accelsearch (Ransom et al. 2002) used to identify po- 
tential pulsars. The original time series were then folded ac- 
cording to the pulse characteristics of each candidate in or- 
der to determine the significance of each signal, construct its 
pulse profile, and refine its characteristics (e.g. period, period 
derivative, etc.) 

We also separately searched the 2008 data for bright in- 
dividual pulses, like those seen from young pulsars such 
as the Crab, using the program single_pulse . py from 
PRESTO, a time-do main matched-filtering techniq ue similar 
to that described by Cordes & McLaughlin (2003). The de- 
dispersed time series were first down-sampled by a factor of 
10 to an effective time resolution of 0.8192 ms. We then con- 
volved the data with square-wave pulses of a variety of widths 
from 1 to 150 samples in duration and searched for peaks sub- 
stantially above those of the noise variations. In no case did 
we find strong pulses that were obviously dispersed in na- 
ture and could be clearly distinguished from RFI. Of course, 
a significantly larger observing bandwidth would provide not 
only more sensitivity to pulses of this sort, but would also dra- 
matically increase the amount of dispersive smearing across 
the band, allowing us to use dispersion to rule out pulses of 
terrestrial origin. While we also searched the 2006 data for 
single-pulse signals without any detections, the sharp jumps 
introduced by the auto-leveler in these datasets rendered it 
particularly difficult to distinguish astronomical pulses from 
system-generated signals. 

3.3. Results 

5 See htt p : / /www, cv. nrao ■ edu/~sransom/presto/| 



Our original search in 2006 yielded candidate pulsed sig- 
nals of high single-trial statistical significance (> 8a) in the 
28-29 June and 10-11 August datasets, independently. The 
highest S/N was obtained at a trial dispersion measure of 
3000pccm~ 3 , and a period of 607 ms in the dataset of 28- 
29 Jun 2006 (Her significance); this pulse profile is shown in 
Fig. |3] Signals at the harmonic periods of 303 ms and 152 ms 
were also weakly detected in this dataset, typically at ~ 4a 
significance. A similar signal, with a period of ~ 604 ms, 
was independently detected in the dataset of 10-11 August, 
with a significance of ~ 8cr at trial dispersion measures of 
3000- 4000 pc cm" 3 ; this is shown in Fig.|4j». Again, har- 
monics of this signal (at 302 ms and 151 ms) were detected at 
lower significance (~ 3.5c). The 607 ms signal was also seen 
on folding the data from 22 June 2006 data using the solu- 
tions determined from the dataset of 28-29 June; this yielded 
a detection significance of 18(7, and the pulse profile of Fig. [5] 
We also tested that these candidate pulsed signals do not arise 
from periodicities in the power jumps in the 2006 data. 

No evidence for a pulsed signal was found in the dataset 
from 2008 at these, or any other, periods. The 2008 data were 
also folded at the pulse periods of the candidates detected in 
2006 (~ 604 -607 ms), without any detectable signal. 

Finally, no other pulsar candidates (above 6a significance) 
were found in any of the datasets. The noise levels achieved 
in each observation are listed in Table Q] The best sensitivity 
was achieved with the datasets of 28-29 June 2006, and 30- 
31 August 2008, which yielded a 10a detection threshold of 
10/iJy, assuming a pulsar duty cycle of 10%. 

4. DISCUSSION AND CONCLUSIONS 
4.1. 77ie 607 ms candidate of 2006 

The high statistical significance of the pulsed signal de- 
tected in the datasets of 2006 implies that it arises either from 
a genuine pulsar towards the GC or as an artifact of a terres- 
trial signal (e.g. RFI). While the non-detection in the 2008 
dataset might indicate the latter possibility, it should be em- 
phasized that the GC environment is very different from the 
environments of typical pulsars. Specifically, a pulsar on a 
short-period (< 100 year) orbit around the GC could easily 
have its emission beam precess away from our line of sight 
over a time-scale of 2 years. This implies that caution must 
be used while dismissing possible pulsar candidates towards 
the GC, although retaining skepticism about their reality. We 
will hence summarize the characteristics of the pulsed signals 
seen in the 2006 datasets, and discuss the possibilities that 
they might arise from a real pulsar or RFI. 

Figs. [3] and [4] show that the detection \ 2 m the two l° n g 
datasets of 2006 increases steadily over the course of each 
observing session, indicating that, if the signals are spurious 
or local RFI, they are at least persistent both over the course 
of each observation, and over multiple observing epochs. It is 
curious, however, that the period of the pulsed signal is dif- 
ferent in the sessions on 28-29 June and 10-11 August; for 
RFI, this would require either that we have detected two dis- 
tinct but alternately intermittent RFI signals, or that the period 
of the RFI itself is changing. Interpreted in terms of Doppler 
shifts, the period change corresponds to a velocity change of 
1500 km/s, much larger than that associated with the Earth's 
motion around the Sun or the motions of terrestrial objects, 
but not implausible for a pulsar orbiting around the GC. 

Next, the S/N ratios of our candidates peak at dispersion 
measures of w 3000-4000pc cm" 3 , comparable to values ex- 
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FIG. 3. — A detection plot of the putative 607 ms pulsar candidate from 
the dataset of 28-29 June 2006; the top panel shows the final pulse profile, 
while the bottom panel shows the growth of reduced \ 2 with increasing data. 
The grayscale denotes the signal strength, with darker signals indicating a 
stronger signal. The white regions correspond to time intervals during which 
the telescope was not pointing at the Galactic Center. 



pected for pulsars at the GC. Unfortunately, our small frac- 
tional bandwidth means that the dispersion in the signal across 
the band is very small (~ 3.6ms over the 800 MHz band- 
width). This means that, unlike the situation in low-frequency 
pulsar surveys, dispersion cannot be used to test whether the 
signal is of extra-terrestrial origin. Finally, the pulse pro- 
files of the candidates are extremely broad, with a duty cy- 
cle of ~ 50%, unlike the narrow profiles expected for high- 
frequency pulsar emission. However, the pulse properties too 
could be affected by the unique GC environment. For exam- 
ple, the thin-screen approximation might not be applicable for 
the scattering, or the screen could be much closer to the pulsar 
than typical estimates of ~ 100 pc; both of these would in- 
crease the scattering time and broaden the pulse profile, even 
at such a high frequency. Specifically, the scattering timescale 
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FIG. 4. — Left: A detection plot of the putative 604 ms candidate from the 
dataset of 10-1 1 August 2006. 

at a frequency of 14.6 GHz for an object at the Galactic Cen- 
ter is 2.5/Z) scat seconds, where D scat is the distance, in pc, of 
the scattering medium from Sgr A*. While the best estimate 
of I>scat is ~ 100 pc from ang ular broadening measurements 
of Sgr A* and nearby masers (Cordes & Lazioiri997f) . the ef- 
fect of scattering material close to Sgr A* is much stronger on 
temporal smearing than on the angular broadening of back- 
ground sources. As such, the angular broadening estimates 
of I>scat do not rule out a substantial contribution to the pulse 
broadening from material closer to Sgr A*. One may hence 
have a sizeable contributi on to the pulse broadenin g from ma- 
terial at D scat < 10 pc (e.g. Macquar t & Bowerl2006l) . The ex- 
pected temporal smearing timescale would then be > 250 ms, 
comparable to that needed to explain the pulse shape of the 
607 ms candidate. The large observed duty cycle of the can- 
didate thus does not rule out the possibility that the signal 
arises from a genuine pulsar. 

It thus appears very difficult to rule out the reality of the 
candidate on the basis of the 2006 data alone, and, as noted 
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FIG. 5. — A detection plot of the pulsed signal on 22 June 2006, after fold- 
ing at the pulse period of the candidate detected in the dataset of 28-29 June 
2006. 



above, the non-detection in 2008 could arise due to preces- 
sion of the pulsar beam away from our sightline. Thus, while 
we remain skeptical about the reality of these signals, we con- 
clude that further observations are needed to test the possibil- 
ity that they arise from a genuine pulsar at the GC. 

4.2. Constraints on the GC pulsar population 

There is compelling but indirect evidence for a substantial 
population of neutron stars at the Galactic Center. However, 
strong interstellar scattering along the line of sight has limited 
past searches for radio pulsars. To overcome these effects, we 
have used the superb sensitivity of the GBT to carry out a deep 
search for pulsars in the central parsec of the GC at 15 GHz 
— the highest observing frequency at which a search has been 
carried out to date. Despite this, we find no convincing pulsar 
candidates. Was our survey sufficiently sensitive to detect a 
population of pulsars around Sgr A*? 

The total number of pulsars detectable at the GC depends on 



the total number of pulsars accumulated in the region, and the 
fraction of these objects that would be detectable given our 
survey sensitivity, and the S/N considerations of Section [2] 
The detectable fraction depends particularly on the number of 
pulsars with flat spectral indices, since these objects influence 
the pulsar luminosity function most strongly at frequencies > 
10 GHz where they are most easily detectable towards Sgr A*. 

A simple estimate of the number of detectable pulsars can 
be obtained by positing that the Sgr A* pulsar population has 
similar properties to those of the known population of pul- 
sars and to estimate the fraction of the known population that 
would be detectable at the GC with our survey. This is done in 
Fig-El where we have plotted pulsars with measu red 1.4 GHz 
luminosities [from the [Manchester et al] d2005l) catalog] on 
a period-luminosity diagram. The solid red line shows the 
pulsar sensitivity curve of our 14.6 GHz survey, obtained us- 
ing equations (1-3) with a 10er detection threshold of 10 /iJy, 
and assuming a 10% pulsar duty cycle, a scattering screen 
distance D scat = 133 pc and a GC dispersion measure 
of 1700pc cm" 3 . The sensitivity curve has been scaled to 
1.4 GHz using a mean pulsar spectral index of -1.7. The cut- 
off in period where most of the sensitivity is lost is taken to 
be at P sp ; n =2 x r scat . This is less severe than the scatter-based 
sensitivity cutoff in Fig.[T]but does reflect the fact that some 
partially recycled or young pulsars (i.e. P< 50 msec) would 
be detectable if they were much brighter than our noise thresh- 
old. 

For comparison purposes, this figure also shows the 
8 .4 GHz sensitivity curve for the Parkes GC survey 
dJohnston et al.i r2006). and the 5cr noise threshold for a deep 
imagin g survey of the GC at 22.5 GHz, using the Very Large 
Array Jzhao et alJ] 2009). again scaling both of these to a fre- 
quency of 1.4 GHz using a mean spectral index of -1.7. We 
also highlight the seven known pulsars with in one degree of 
Sgr A*, including four new ones from iDeneva et al.1 (2009) 
and lCamilo et alJ 00091) . 

A more rigorous estimate can be obtained by comput- 
ing the fraction of pulsars detectable above some flux den- 
sity cutoff by considering the pulsar luminosity function 
at vo = 1.4 GHz, fo(L), combined with the spectral in- 
dex distributi on, p( a). This is the app roach followed by 
iPfahl&Loebl (120041) and lCordes & Laziol (119971) . but updated 
with the most recent results on pulsar luminosity fun ctions 
dFaucher-Giguere & Kaspil 120061: lLorimer et ail |2006j). We 
model the 1.4 GHz luminosity function as a power law be- 
tween lower and upper cutoffs L m ; n and L max respectively: 



fo(L)=AL- 



A = (1- ( S) 



- 1-/3 _ J l ~P 



(3) 



where the normalization is chosen so that the integral over 
all luminosities is unity, such that fdL is interpreted as 
the fraction of all pulsars with luminosities between L and 
L + dL. Recent studies suggest t hat L m j n = 0.01 mJykpc 2 , 
L ma i = 32 Jy kpc 2 and = 1.2-2 dFaucher-Giguere & Kaspil 
2006; Lorimer et al]l2006l) . Following|SmitsetalJ 42009J), the 
spectral index distribution is modeled as a Gaussian 



p(a) = 



1 



exp 



(a-a m y 



2al 



(4) 



with mean spectral index a = -1.7 and standard deviation 
a a = 0.35. 

The pulsar luminosity function at some arbitrary frequency 
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FIG. 6. — The 1.4 GHz luminosities of the known sample of pulsars versus pulse periods (blue dots). Larger circles (red dots) indicate those pulsars within one 
degree radius of Sgr A*. The IOct pulsar sensitivity of our 14.6 GHz search is shown by the solid red line. This was obtained by using the flux density limit (10 
/^Jy at IOct) of our survey to calculate the luminosity limit at the distance of the Galactic Center and then scaling the result to a frequency of 1.4 GHz, using an 
average spectral index of (a) =-1.7. We also sh ow sensitivity curve s derived in the same manner for a Parkes 8.4 GHz pulsar survey of the GC (the dashed red 
line, with a 10<r dete ction threshold of 200 /ijy; [Johnston et al. 2006), and a deep 22.5 GHz VLA image of the GC (the dotted green line, with a 10ct detection 
threshold of 200 ^Jy; Zhao et al. 2009). 



is then 



f(v,L)= dap(a)f[L 



(5) 



where, for a given spectral index a chosen from the distribu- 
tion, the luminosity function / has lower and upper cutoffs 



V ■ -T 

^min — Mmn 



V 



V -T 



V 



v j \ v a 

and the normalization constant A is modified to 



(6) 



(7) 



in order to ensure that fdL may be interpreted as the fraction 
of all pulsars in the luminosity range L to L+dL. We inte- 
grate equation (O to obtain the total fraction of pulsars above 
a given flux density threshold L cut = d 2 S C ut- Of course, this 
threshold depends on both the pulsar period and observing 
frequency due to propagation effects and changes in the sys- 
tem temperature. Figure [7] hence plots the fraction of the pul- 
sar population detectable above 10<7 significance at 15 GHz 
versus pulse period. 

Note that the present survey is sensitive only to the slower 
(P> 40 ms) and more-luminous (L>40 mJy kpc 2 ) pulsars. 
The millisecond pulsars that are pres umably poweri ng the 
low-mass X-ray binaries near Sgr A* (Muno et al. 2005) and 
the low-luminosity tail of young pulsars (Camilo et alj |2009t) 
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FIG. 7. — The fraction of Galactic Center pulsars that would be detectable at 10cr significance for a 10 h integration on the GBT at 15 GHz as a function of spin 
period. We assume an intrinsic pulse width of 10% of the spin period, and take D 8ca t = lOOpc. The curves, from top to bottom, correspond to luminosity function 
indices of = 1.2, 1.5 and 1.8, with a mean spectral index of a = —1.6 and dispersion <r Q = 0.35, as discussed in the text. 



are out of the reach of our survey. However, it is also clear 
from the figures that this is the first survey capable of peer- 
ing past the "fog" of scattering material and detecting a sig- 
nificant number of pulsars within a parsec of the GC with 
properties similar to the kn own pulsar population . Past high- 
frequency pulsm_S£arches (Johnst on et alj 12006) or imaging 
searches dZhao et al.H2.Q09h have not had the requisite tempo- 
ral or flux density sensitivity to detect a significant fraction 
of the known population. The GBT search thus represents a 
significant improvement over past pulse searches and imag- 
ing efforts. We note, in passing, that this implicitly assumes 
that pulsar spectral indices do not typically steepen at high 
frequencies, > 5 GHz. 

It is clear from Figs.[6]|7]that the 15 GHz GBT survey could 
have detected a significant fraction (~ 1%-15%) of the pul- 
sars around Sgr A*, if they had properties similar to those of 
the known population. The estimate obtained from Fig. [6] is 
at the high end (15%) and it is likely biased by luminosity- 
dependent completeness limits in pulsar surveys. The low- 
est estimate (~1%) comes from the curve in Figs.|7]with the 
steepest luminosity slope (/? = 1 .8). We adopt a nominal value 
of 5% from the intermediate curve ((3=1 .5) which is based on 
our best current knowledge of the properties of the pulsar pop- 
ulation and the scattering material toward Sgr A*. Given this 
detection fraction and our null detection we can use straight- 
forward binomial statistics to estimate the size of the putative 
pulsar population at the GC. If the probability of detecting a 
normal pulsar is 5%, the non-detection of any pulsars in our 
survey implies the upper limit (at 99% confidence level) of 
90 normal pulsars within the 1 pc region around the GC en- 
compassed by the GBT beam. Taking the full allowed range 
of the detection fraction (1% to 15%), the upper limit on the 



number of normal pulsars ranges from 460 to 30, respectively, 
again at 99% confidence level. Although our estimate is both 
approximate and subject to much uncertainty, we note that 
it i s significantly lower than the ~ 100- 1000 pulsars derived 
by dPfahl & Loebl2004l) for the normal pulsar population with 
orbits of < 100 yrs (i.e. a radius 50 times smaller than the size 
of the GBT search area). 

Finally, we have shown that the frequency range 10- 
16 GHz is optimal for searches for "normal" pulsars at the 
GC. The GBT remains the most powerful high-frequency in- 
strument capable of detecting the GC pulsar population for 
at least the next decade, until the advent of next-generation 
telescopes like the Square Kilometer Array. The primary lim- 
itation of the GBT (and of the present survey) is the relatively 
small instantaneous bandwidth (800 MHz) available for such 
searches, resulting in a small fractional bandwidth. An in- 
creased bandwidth at the GBT would imply not only an im- 
provement in sensitivity, but also a better rejection of terres- 
trial signals, using the dispersive sweep of genuine signals 
across the band. Fig. [6] shows that an improvement in sensi- 
tivity by merely a factor of 2 - 3 would push the GBT into 
the bulk of the pulsar population. Future GBT experiments 
should hence aim to utilize the full frequency coverage avail- 
able with the high-frequency receivers, for both better dis- 
crimination against systematic effects and improved sensitiv- 
ity. 
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